Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 14 August 2012 



(MN IMfeX style file v2.2) 



A Refined Measurement of the Mean Transmitted Flux in 
the Lya Forest over 2 < z < 5 Using Composite Quasar 
Spectra 



George D. Becker 1 *, Paul C. Hewett 1 , Gabor Worseck 2 and J. Xavier Prochaska 2 

1 Kavli Institute for Cosmology and Institute of Astronomy, Madingley Rd, Cambridge, CB3 0HA, UK 

2 Department of Astronomy and Astrophysics, UCO/Lick Observatory, University of California, 1156 High Street, Santa Cruz, 
CA 95064, USA 



Draft version 14 August 2012 



Of 



ABSTRACT 

We present new measurements of the mean transmitted flux in the Lya forest 
over 2 < z < 5 made using 6065 quasar spectra from the Sloan Digital Sky Survey 
DR7. We exploit the lack of evolution of the mean quasar continuum to avoid the 
bias introduced by continuum fitting at high rcdshifts, which has been the primary 
systematic uncertainty in previous measurements of the mean Lya transmission. The 
individual spectra are first combined into twenty-six composites with mean redshifts 
spanning 2.25 ^ z comp ^ 5.08. The flux ratios of separate composites at the same rest 
wavelength are then used, without continuum fitting, to infer the mean transmitted 
flux, F(z), as a fraction of its value at z ~ 2. Absolute values for F(z) are found by 
scaling our relative values to literature measurements made from high-resolution data 
at z ^ 2.5, where continuum uncertainties are minimal. We find that F(z) evolves 
smoothly with redshift, with no evidence of a previously reported feature at z ~ 3.2. 
This trend is consistent with a gradual evolution of the ionization and thermal state 
of the intergalactic medium over 2 < z < 5. Our results generally agree with the most 
careful measurements to date made from high-resolution data, but offer much greater 
precision and extend to higher redshifts. This work also improves upon previous efforts 
using SDSS spectra, which have largely suffered from systematic errors. This is the 
most accurate and precise measurement of the mean transmitted flux at these redshifts 
to date, and should serve as a reference for future studies that use the Lya forest. 

Key words: intergalactic medium - quasars: absorption lines - cosmology: observa- 
tions - cosmology: large-scale structure of the Universe 



1 INTRODUCTION 

The pattern of absorption lines imprinted on the spectra 
of distant objects by neutral hydrogen in the intergalactic 
medium (IGM), known as the "Lya forest", is one of the 
most fundamental probes of cosmic structure. The charac- 
teristics of this absorption reflect the density distribution, 
ionization state, and temperature of the intergalactic gas. 
The gas, in turn, closely traces the underlying distribu- 
tion of dark matter, albeit with significant deviations arising 
from hydrodynamics coupled to the radiative and mechani- 
cal feedback from galaxies and AGN. Consequently, the Lya 
forest allows us to probe the baryon physics connecting the 
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evolution of large-scale structure to the highly non-linear 
processes driving galaxy formation. 

A key requirement for utilizing the Lya forest to study 
either cosmology or baryon physics is to accurately es- 
tablish the basic parameters describing the flux distribu- 
tion. The most basic of these is the mean transmitted flux, 
F = (e~ T ), and its evolution with redshift. The mean flux, 
which is often expressed in terms of an effective optical 
depth, T c ff(z) = — hxF{z), directly constrains, for exam- 
ple, the intensity of the metagalactic ionizing background 
(e.g-lRauch et al.lll997l : iMcDonald fc Miralda-EscudeboOll ; 
iBolton et alj 12003 ). It also influences statistics such as the 
flux probability distribution function and power spectrum, 
among others, which are used to constrain qu antities rang 
ing from the temperature o f the IGM (e.g., 
l2000l : IZaldarriaga et al.ll200ll : iLidz et al.|[201ol : 
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2011) to the free-stre aming length of dark matter particles 
fe.g.. IViel et al.l [20081 ). Lya forest studies increasingly rely 
on comparing real data to artificial absorption spectra drawn 
from simulations, which must be calibrated to the correct 
mean flux. With the increasing availability of high-quality 
quasar and gamma ray burst spectra taken with large tele- 
scopes, and the large number of moderate-resolution quasar 
spectra gather ed with surveys s uch as the Sloan Digital Sky 
Survey CSDSS; lYork et al.ll200C) and the Baryon Os cillation 
Spectroscopic Survey (BOSS; Dawson et al J 120121) . a pre- 
cise measurement of F(z) is therefore critical for accurately 
extracting the maximum amount of information from the 
forest. 

The main challenge in measuring F(z) is to overcome 
cosmic variance while accurately estimating the unabsorbed 
continuum of the background objects, which have so far 
typically been quasars. Studies using high-resolution (gen- 
erally R ~ 40, 000) spectra have the advantage that the 
peaks of the Ly a forest transmission should approach the 



continuum (e.g., Rauch et al.l 1 19971: IMcDonald et al] 



Schave et all [2003: 



ISong aila 



2004: 



iKirkman et al.l 



200 



200 



Kim et al.l 120071; iBecker et al. 120071; IFaucher-Giguere et al.l 



20081 : iDall'Agho et all l20Qg| ; iBecker et all 1201 If ). Even at 
moderate redshifts [z ~ 2 — 3), however, voids in the IGM 
will have a non- negligible optical depth. The continuum will 
therefore lie somewhat above the tops of the transmission 
peaks, a bias that increases towards higher redshifts. Al- 
ternatively, studies using large sets of moderate-resolution 
spectra have either attempted to simultaneously solve for 
F(z) and the quasar continua in a strictly statistical sense 
|Bernardi et al.ll2003l ). or have fit continua on an individual 
basis using principle co mponent analysis |McDonald et all 
120051 : iParis et alj|201l|j_ or with bias corrections estimated 
from model spectra jDaH'Aglio et all 12009] ). At the high- 
est redshifts (z > 5), F(z) measurements have generally re- 
lied on extrapolating a po wer-law continu u m from redward 
of the Lya emiss ion line jSongailal |2004 iFan et al] l200rj ; 
IBecker et ai]|2007h . 

The most accura te F(z) measurement t o date has ar- 
guably been made by IFaucher-Giguere et all (|2008l ). These 
authors used a set of 86 high-resolution quasar spectra and 
made statistical corrections to the continua as a function 
of redshift based on artificial spectra drawn from hydrody- 
namic simulations. Notably, they identified a possible sharp 
(Az < 0.5) dip in the evolution of the me an opacity near 
z ~ 3.2. A similar feature was found by iBernardi et al] 
(|2003l ) using a very different data set and measurement 
technique. The optical depth of the optically thin IGM 
to Lya will scale with the H I neutral fraction, which in 
turn depends on the gas temperature, T, and the H I ion- 
ization rate, F, as /hi oc y-o^p- 1 ^ where the temper- 
ature dependence is for case A recombination. Bernardi 
et al. noted that a decrease in t c h(z) at z ~ 3.2 could 
indicate a temporary increase in the IGM gas tempera- 
ture accompanying the end of helium reionizatio n (see also 
iTheuns et aT]|2002l : IFaucher-Giguere et al]|2008l ). However, 
the reality of this feature remains unclear. Although a simi- 
lar dip was found at modest significance bv lDalFAglio et ail 
(2008 1) m a sample of 40 high-resolution spectra , neither 
IMcDonald et al] (|2005l ) nor lDaU'Aglio et al] j2009l ) detected 
such a feature using sampl es of moderate-resolut ion SDSS 
spectra larger than that of IBernardi et all (|2003l ). Indeed, 



IFaucher-Giguere et al] (|2008l ) emphasized that the signifi- 
cance of their feature may have been exaggerated by the 
assumption of an underlying power law evolution in r e ff. 
The association of such a sharp feature with helium reion- 
ization has a lso b een challenged on theoretical grounds. 
iBolton et all l|2009l ) argued that even if the IGM is photo- 
heated during the course of a very rapid (Az ~ 0.2) helium 
reionization, the expected signature would be a sudden de- 
crease in T e s followed by a gradual recovery driven by adi- 
abatic cooling in the voids. Moreover, recent temperature 
measurements indicate that heating due to hel ium reioniza- 
tion was an extended process with Az > 1 (|Becker et al] 
l201lfh and so it is unclear why such a sharp feature in the 
opacity evolution should occur. 

In this paper we present a new and highly precise mea- 
surement of the mean transmitted Lya flux over 2 < z < 5 
using a sample of 6065 moderate-resolution quasar spectra 
drawn from the seventh data release of the Sloan Digital 
Sky Survey (SDSS DR7). The main innovation of this work 
is the introduction of a new method that uses composite 
spectra to measure the overall shape the mean flux evolu- 
tion to high accuracy without fitting continua. These re- 
sults are then normalized to measurements made from high- 
resolution data at z ~ 2, where continuum errors are mini- 
mal. Our approach builds on the successful use of composite 
quasar spectra to statistically constrain other properties of 
the I GM, for example the mean free path of ionizing pho- 
tons (|Prochaska et al]|2009l ). We also take advantage of the 
improved statistical accuracy afforded by DR7 to increase 
substantially th e number of quasars analyzed compared to 
previous works (IBernardi et alj|2003l; IMcDonald et al.|[2005l : 
IDairAglio et al]|2009l : |Paris et al]|201ll . although Paris et al. 
also used DR7 data). This paper should therefore serve as a 
reference for future studies that utilize the mean Lya opac- 
ity. In Section [2] we describe our technique, which uses com- 
posite spectra to compute F(z) relative to that at z ~ 2. Our 
results are presented in Section [3] In Section |4j we compare 
our measurements to results from the literature. Finally, we 
summarize in Section [S] 



2 DATA & METHOD 

The normalized transmitted flux in a single spectrum is 
given by the observed flux divided by the unabsorbed con- 
tinuum, F 1 — /oba//cont- The goal of this paper is to mea- 
sure the mean normalized flux as a function of observed 
wavelength (i.e., redshift) across our entire sample of ob- 
jects, F(z) — (F'(z)}, where 1 + z = A b s /ALya, and 
Xhya = 1215.67 A. Rather than fit continua to individual 
spectra, however, we use the fact that the mean continuum 
for a large group of quasars will be highly similar regardless 
of redshift. The difference between composites at different 
redshifts should therefore primarily reflect the evolution in 
the mean transmitted flux in the Lya forest. In this section 
we describe how our composites are constructed, and how 
these are used to calculate F(z). 



2.1 Constructing composite spectra 

We constructed 26 composite quasar spectra using a total 
of 6065 individual spectra from the SDSS DR7 quasar cata- 
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Composite spectra used in this work. Columns give the 
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4.346 


70 
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4.495 


101 
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1.8 


4.694 


63 


1.8 


5.0 


4.896 


51 
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5.2 


5.083 


36 



logue |Schneider et al.ll2010l ). The redshift range, mean red- 
shift, and the number of quasars contributing to each com- 
posite are listed in Table [T] Redshift bins of Az = 0.1 were 
used between z = 2.2 and 4.4. Above z = 4.4 we used a bin 
size of Az = 0.2 in order to increase the number of objects in 
each composite. The quasars were selected to have an abso- 
lute magnitude Mi > —29.0, while lacking broad absorption 
line (BAL) features. For quasars present in the SDSS DR6 
rel ease, BAL quasars were defined using the classifications 
of lAllen et alj (|201lf ). while the small fraction of quasars 
present only in DR7 were classified via visual inspection. 
None of the results presented here, however, are sensitive 
to the exact definition of the BAL quasar sample. We also 
excluded objects where the peak of the Lya emission line 
was greater than ten times the continuum level near a rest 
wavelength of 1280 A. This was done primarily to avoid very 
high-redshift objects with extremely low signal-to-noise ra- 
tios in the continuum (S/ N < 2). Improved quasa r redshiftsQ 
were determined using the lHewett fc Wild! (|2010l ) scheme for 
generating self-consistent redshifts for quasars with SDSS 
spectra. We used spectra that had b een corrected for tel- 
luric emission line residuals in the red (|Wild fc Hewettll2005l. 
2010). In addition, for roughly 6 per cent of our objects, 
spectra from multiple epochs were combined to increase the 
signal-to-noise ratio. 

The SDSS spectra are supplied on a logarithmic wave- 
length scale, with pixels of 69 km s _1 width. When generat- 
ing composite spectra in the quasar rest-frame, nearest-pixel 
values were employed without interpolation or rebinning. 

1 http:/ /www.sdss.org/dr7/products/value_added/#quasars 
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Figure 1. Continuum power law slopes for the quasars used to 
construct our composites, where fx <x X a , and a is measured 
between 1270 A and 1510 A (rest wavelengths). Black points give 
the values for individual objects. Red points with error bars show 
the mean value and standard deviation in each redshift bin. The 
horizontal line at a = — 1.56 shows the composite value from 
IVanden Berk et al.1 l|200ll V 



The arithmetic mean over all contributing objects was then 
computed at each pixel. All spectra were weighted equally 
in order to minimize the effects of cosmic variance on the 
flux in the Lya forest. Error arrays were also computed by 
directly computing the intra-sample standard deviation. For 
the initial composites (described below), pixels falling within 
16 A intervals centered on the wavelengths of strong [O i] sky 
emission lines at 5578.5 and 6301.7 A were masked out. For 
the final composites, we masked all regions contaminated 
by strong telluric lines that are flagged in the bitmask ac- 
companying each SDSS spectrum. We note that the flagged 
regions typically do not included all of the pixels in a spec- 
trum affected by skylines, and often cover only the core re- 
gion contaminated by a skyline and not the line wings, which 
can still be significantly affected by residuals. Using spectra 
that have been corrected for skyline residuals therefore still 
proved beneficial. 

Initial composites were generated using the original 
quasar spectra without modification, except to divide each 
object by its mean flux near 1285 A, which is the region of 
the continuum closest to the Lya forest that is relatively 
free of emission lines. These composites exhibited an en- 
couraging similarity in both continuum and emission-line 
properties longward of 1250 A in the quasar rest-frame. A 
weak but systematic trend, however, appeared in the over- 
all spectral shape as a function of redshift. To correct for 
this, we re-calculated the composites after adjusting the con- 
tinuum slope of each individual object using the following 
procedure. The high signal-to-noise ratio initial composite 
at z ~ 2.65 was chosen as a reference. For each individual 
quasar, a power law was then fit to the ratio of the flux in its 
spectrum to that in the reference composites, and the spec- 
trum was divided by this fit at all wavelengths. The fit was 
made between rest-frame wavelengths 1270 A and 1510 A, 
which is the reddest rest-frame wavelength region covered 
in the SDSS spectra for objects at z ~ 5. Note that we 
avoid the high-ionization lines Lya A1216, N v A1239,1243, 
and C IV A1548,1551, which tend to vary strongly between 
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Figure 2. The twenty-six composite quasar spectra used in this work, plotted verses rest-frame wavelengths. The spectra are binned 
using 345 km s -1 pixels, as used in the analysis, and separated into multiple panels for clarity. The mean quasar redshift is also displayed 
for each composite. The vertical dotted lines indicate the region of the Lya forest used to measure F(z). Note that the prominence of 
the emission lines in the forest naturally decreases with increasing composite redshift as the overall transmitted flux decreases. 



objects. We also exclude the region between 1380 A and 
1410 A, which contains the Si iv A1394.1403 emission line 
doublet. 

The best-fitting spectral indices, calculated by taking 
the slope of the initial z ~ 2.65 composite as a reference 
point, are plotted in Figure [T] We note that higher-redshift 
objects tend to be systematically redder, albeit slightly. This 
trend may reflect a combination of effects related both to the 
selection of the quasars in the SDSS and intrinsic changes 
in the quasar spectra, or even aspects of the data reduction. 
It is clearly essential, however, to remove any systematic 
trend that could impact the determination of the measure- 
ment of the Lycv-forest depression with redshift. The only 
assumption underlying the nature of the correction proce- 
dure adopted is that the essentially featureless changes to 
the continuum slopes, which were measured over large wave- 
length intervals redward of the Lya emission, extend into the 
spectral region used to study the Lya forest. We test this 
assumption in Section [3] 

The final composites generated from the continuum- 
corrected SDSS spectra are shown in Figure [2] The compos- 
ite spectra are strikingly similar redward of the Lya+N v 
emission lines. The overall amplitudes and slopes are ex- 
pected to be the same based on the continuum correction 
procedure described above. We note, however, that the am- 
plitudes of the emission lines redward of Lya+N v are also 
very similar. This is particularly true of the low ionization 



lines, O I A1302 + Si n A1304 and C n A1334, which is 
encouraging for our purposes since the lines in the forest 
are a lso from low ionization species (e.g. IVanden Berk et al.l 
l200lD . 

Significant variations do appear in the strength of the 
Lya emission line, which, taken at face value, might be of 
concern. The effect, however, is due to errors in the redshifts 
of the individual objects, which, because of Lya absorption 
on the blue side of the Lya emission line, tend to decrease 
the height of the emission in the composite spectrum. The 
effect increases towards higher redshift as the absorption in 
the Lya forest increases. The redshift uncertainties for the 
S DSS quasar spectra al so increase with redshift (Table 3 
of iHewett fc Wild! 1201(f ) as the C IV emission line has to 
be used in the redshift determination. At the very highest 
redshift, 2 comp ~ 5.08, the effect apparently diminishes (see 
Figure [3]) because the Lya emission line itself dominates in 
the redshift determination, causing the peak of the compos- 
ite Lya line to sharpen by selection. By comparison, the flux 
over the region spanning 1250 A < A < 1500 A, which con- 
tains relatively weak emission lines, shows very little varia- 
tion with redshift. Since the the emission lines in the forest 
are also weak, redshift errors should have a minimal effect 
there, as well. 

We directly verified that redshift uncertainties affect 
the strength of Lya by comparing composite spectra of 56 
quasars at 4.4 < 2 < 5.2 with improved redshifts from 
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Figure 3. A demonstration of the effects of rcdshift errors on 
composite spectra. The same set of 56 4.4 < z < 5.2 quasar spec- 
tra were used to create both versions of the composite plotted 
here. The red line shows the composite as generated using in- 
dividual rcdshifts measured from the SDSS spectra themselves, 
while for the black line the redshifts were measured from higher- 
quality GMOS spectra. The improved rcdshifts from the GMOS 
data produce a more sharply peaked Lya emission line but have 
little effect on the rest of the spectrum. Note that, since we se- 
lected only objects where the GMOS redshifts could be measured 
from O I A1302 + Si II A1304 and/or C II Af334, the emission lines 
are somewhat stronger than the composites shown in Figure [2] 



recently-obtained Gemini/GMOS spectra (Worseck et al, in 
prep). Two composites were generated using the same set 
of SDSS spectra. In one case we use the redshifts measured 
directly from the SDSS spectra, while in the other case we 
use redshifts from O I A1302 + Si II A1304 and/or C n A1334 
measured from the GMOS data. The results are shown in 
Figure [3] The improved GMOS redshifts produce a stronger 
and more sharply peaked Lya line, as expected, while the re- 
mainder of the spectrum is nearly unaffected. This supports 
our assumption that, for all of our composites, the differ- 
ences in the fluxes blueward of the Lya emission lines are 
due solely to the evolution of the mean opacity of the Lya 
forest. It is these differences that we will use to constrain 
F{z). 

2.2 Calculating the mean transmitted flux 

At a given rest wavelength Ao , the flux of a quasar composite 
with mean emission redshift Zcomp will be attenuated by a 
factor 

F(z) = /° b °( A °-^°'°p) | (1) 

/cont(Ao) 

where 1 + z = Ao(l + z conl p)/ Xhya- We unfortunately do not 
have a completely unabsorbed composite at 2 comp = that 
would allow us to measure / CO nt(Ao) absolutely, and hence 
calculate F(z) directly. However, we can measure the ratio 
of F(z) at two redshifts by evaluating the ratio of fluxes 
between two composites at the same rest-frame wavelengths, 

F{Z\) _ /obs(Ao, Zcomp.l) ^) 
F{Z2) /obs(A(), Zcomp,2) 



Here, the dependence on the continuum has been eliminated 
because we are using composites rather than individual spec- 
tra. This relation will hold to the degree that the mean con- 
tinuum in the composites either does not evolve or has been 
adequately corrected with redshift, an assumption we test in 
Section^ Nominally, however, this dataset allows us to com- 
pute F(z) up to a normalization factor by finding a function 
that satisfies the set of observed flux ratios. 

We parametrized F(z) as a stepwise function in bins 
of Az = 0.1, with the first bin centered on z — 2.15. We 
then computed F(z) as a fraction of the transmitted flux in 
that bin, i.e., F(z)/F(z = 2.15). When analyzing the data, 
we bin the composite fluxes into 345 km s _1 pixels. This is 
to avoid large non-symmetric errors in F{z\) / F{z2) when 
F(z2) becomes small, although nearly identical results were 
obtained without binning. We also made secondary adjust- 
ments to the composite continua by dividing each composite 
by a power-law fit over 1270 A < A < 1510 A to the ra- 
tio of that composite to the average of all composites with 
2.25 < Zcomp < 4.25. This was a minor correction to help 
remove any residual differences in the composite shapes, and 
changed the final F(z) values by less than the la errors. 

The primary goal of this paper is to measure the mean 
opacity of intergalactic hydrogen to Lya. Here we loosely 
define "intergalactic" to include all diffuse, highly-ionized 
absorbers with H I column densities A^hi < 10 19 cm -2 . As 
defined in equation [TJ however, F(z) in the Lya forest will 
also include absorption from metal lines, as well as from 
higher-column density Lya absorbers. We now describe how 
we account for these additional contributions to arrive at 
the desired F(z) values. 

2.2.1 Metal lines 

Absorption in the Lya forest due to intervening metals, 
while small compared to the Lya absorption, must be ac- 
counted for in order to measure F(z) precisely. Conventional 
studies using cont inuum- normali zed spectra hav e relied on 
mask ing (e.g.. ISchave et al"1l2003t ) or fitting out (Ki m et al] 
l2007t ) metal lines indi vidually, or else apply ing global statis- 
tical corrections (e.g.. iKirkman et al]|2005T ). Measuring flux 
ratios in composite spectra, however, offers an alternative 
approach to mitigating the impact of metals, as discussed 
below. 

The mean level of metal-line absorption may vary, in 
principle, with both quasar redshift and rest-frame wave- 
length. A large fraction of absorption in the forest will 
come from the C IV A1548, 1551, Si IV A1394, 1403, and 
Mg n A2796, 2804 doublets, which, because they have rest 
wavelengths longward of Lya, will produce, on average, a 
continuum of absorption extending over the region redward 
of the Lya emission line we use to scale our spectra and 
into the forest. An additional contribution will come from 
the various metal lines associated with high H I column den- 
sity systems such as DLAs. These lines have a variety of rest 
wavelengths, including some that fall between Ly/3 and Lya. 

In practice, the flux decrement from metals is relatively 
minor, and appears to depend on ly mildly on e ither quasar 
redshift or rest-frame wavelength. iBecker et al] (|201lf ) found 
a nearly constant decrement of ~3 per cent in the region 
between the N v A1239, 1243 and Si IV A1394, 1403 emis- 
sion lines in the high-resolution spectra of a large sample 
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of quasars with redshifts z ~ 2 — 4. iKirkman et alj (|2005l ) 
also found that the absorption due to metals redward of 
Lya evolves s lowly with quasar redshift. In the Lya forest, 
ISchave et all (|2003h found a similar (~3 percent) level of 
absorption from metals over Lya redshifts z ~ 2 — 4, with 
little evidence for strong evolution. iKim et al.l (|2007l ). in pos- 
sibly the most careful such work to date, also found a mean 
decrement of ~3 per cent due to all metals in the forest over 
2 ~ 2 — 3, albeit with significant scatter between sight lines, 
and no evidence for evolution. 

In a composite spectrum, the absorption due to met- 
als will essentially scale the mean flux by a multiplicative 
factor, and may also introduce a change in slope with rest- 
frame wavelength. Given the apparently mild dependence of 
this absorption on either rest-frame wavelength or quasar 
redshift, however, we expect these effects to be minor (on 
the order of a few percent). More importantly, the process 
of scaling the composites to have the same amplitude and 
slope between the Lya and C IV emission lines will tend 
to eliminate the relative differences in the metal line ab- 
sorption as a function of composite redshift. This is another 
advantage of our differential approach to measuring F(z). 
Systematic errors with redshift may remain if the metal- 
line absorption (or more precisely, the relative difference in 
metal-line absorption between composites at different red- 
shifts) is not well represented by a power law in wavelength 
(see Section |2.1|) . Given the metal-line measurements dis- 
cussed above, however, we expect such residual systematics 
to be minor, even compared to relatively small statistical 
errors on F(z) from such a large data set. 

We note that establishing an absolute scaling for F(z) 
still requires a correction for metal lines, at least over a 
narrow range in redshift. As described below, we use abso- 
lute values_jbr_Jj , (2_~ 2) that have been corrected follow- 
ing [Scha^eFlil] |2003h . This correction for metal lines is 
roughly consistent with other works, although it is still a 
potential source of systematic error in our final F(z) values, 
albeit at the <1 per cent level. 

2.2.2 Correction for optically thick absorbers 

The total Lya opacity will include a contribution from ab- 
sorbers that are optically thick to ionizing photons (i.e., with 
Nhi ^ 10 2 cm -2 ). As discussed below, most of this ab- 
sorption will come from super Lyman limit systems (SLLSs; 
10 190 cm~ 2 sC TVhi < 10 20 ' 3 cm" 2 ) and damped Lya ab- 
sorbers (DLAs; Nhi ^ 10 20,3 cm~ 2 ). These are believed to 
be associated with intervening galaxies, and are difficult to 
model in numerical simulations. In order to facilitate a direct 
comparison with simulations, therefore, we correct for these 
abs orbers in order to arrive at ou r final F(z) values. Follow- 
ing [^Mcher^Giguire^r|al] (|2008l ). we use a model to calcu- 
late the integrated absorption due to optically thick systems. 
For a given frequency distribution in column density, Nhi, 
Doppler parameter, b, and redshift, the flux decrement due 
to these absorbers will be 

l-JBr(a) = 1±£ r^dNm f db f{N m ,b, z)W {N m ,b) ,(3) 

where Wo is the rest-frame equivalent width. At high col- 
umn densities, Wo will depend mainl y on the column density, 
and so we set b = 20 km s" 1 , as in IFaucher-Giguere et all 



|2008l ). We assume that the shape of the column den- 
sity distribution does not evolve with redshift, such that 
f(Nm,z) = f (Nhi) dn/dz. We use the broken power- 
law column density d istribution at z ~ 3.7 compiled by 
IProchaska et alj (|2010h . and adopt 10 19 cm" 2 and 10 22 cm" 2 
as our limits of integration. Lyman limit systems with 
10 17 ' 2 cm" 2 < TVhi < 10 19 cm" 2 were found to contribute 
a negligible amount of Lya absorption at our level of pre- 
cision, while DLAs with Nhi > 10 22 cm" 2 are vanishingly 
rare. The redshift evolution of the number density of SLLSs 
and DLAs is not known precisely over the entire redshift 
range; however, dn/dz oc (1 + z) 2 gives a reasonable fit 
to the evolution of DLAs o ver 2.3 < z < 4.4 measured 
bv lProchaska fc Wolfel l|2009h . Adopting this evolution gives 
an integrated flux decrement due to SLLSs and DLAs of 
l-F T (z) = 0.0045 [(l + z)/3] 3 '°. We note that this is some- 
what more absorption than estimated by Faucher-Giguere 
et al., who adopted a model with f(Nui) oc N-^l' 5 over all 
column densities, and dn/dz oc (1 + z) . 

When calculating F(z), we first fit raw values relative 
to z ~ 2 that include absorption from both the optically 
thin Lya forest and from optically thick absorbers. We then 
corrected these relative values using the above scaling for 
Ft to give the relative mean flux in the forest only. These 
corrected values of F(z)/F(z = 2.15) were then scaled to 
measurements at z 2.5 from the literature to derive abso- 
lute values for F(z), as described below. 

2.2.3 Absolute scaling 

In order to compute absolute values for F(z), we scaled 
our results to F(z) measurements made from high-resolution 
spectra at redshifts where co ntinuum uncertainties are min - 
imal. We used the results of IFaucher-Giguere et all (|200Sf ) . 
who used artificial spectra drawn from cosmological simula- 
tions to correct for continuum-fitting biases. These authors 
found the corrections to be small (< 2 per cent) at z < 2.5 . 
Their results were also confirmed by iBecker et al.l (|201lf ) , 
who measured F(z) from high-resolution spectra using a dif- 
ferent approach to correct for continuum errors. Becker et 
al. re-normalized both the real data and the simulated spec- 
tra by deliberately placing the continuum along the tops of 
the highest transmission peaks. The true F(z) was then mea- 
sured from the simulations once the simulated optical depths 
had been adjusted such that the re-normalized data and sim- 
ulations had the same mean flux. The resulting F(z) agreed 
with the Faucher-Giguere et al. measurements to within one 
per cent over 2 < z < 2.5. We note that Faucher-Giguere 
et al. and Becker et al. also used different methods to cor- 
rect for metal-line contamination; the Faucher-Giguere et 
al. values adopted here include a corre ction based on metal- 
line measurements in the forest from ISchave et all (|2003h . 
whereas Becker et al. applied a similar correction based on 
absorption redward of Lya. 

To scale our measurements, we interpolated our 
F(z)/F(z — 2.15) results onto the centers of the four 
Az = 0.1 redshift bins betwee n z = 2.2 and 2.5 used 
by IFaucher-Giguere et all (120081 1. The mean flux ratio be- 
tween our raw results and the Faucher-Giguere et al. values, 
weighted by the errors, was then used to compute our final 
F(z). We note that when computing the scaling we first ad- 
justed the Faucher-Giguere et al. data to account for the dif- 
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ference in the amount of absorption we assume for optically 
thick absorbers. This is a minor (~0.1 per cent) correction 
at z ^ 2.5, however. 

2.2-4 Error estimates 

Uncertainties in F(z) were estimated using a bootstrap re- 
sampling approach. In each iteration, a new composite was 
generated for each redshift bin from a set of Af quasars 
drawn randomly, with replacement, from the set of Af ob- 
jects in that bin. F(z) was then re-computed from the new 
composites following the procedure described in Section [2T] 
This process was repeated 1000 times, and the full covari- 
ance arrays for both the unsealed F(z)/F(z = 2.15) and 
scaled F(z) were computed from the results. We stress that, 
because we are fundamentally fitting ratios of fluxes at dif- 
ference redshifts, the errors in F(z) will be correlated. Any 
analysis that uses the F(z) results should therefore consider 
the entire covariance array. 

2.2.5 Impact of sky-subtraction errors 

Measuring F(z) at z < 2.3 requires analyzing fluxes at very 
blue wavele ngths (Aobs < 4000 A). Recent papers using 
SDSS data (|Prochaska et alj|2010l ; iParis et al]|201ll ). how- 
ever, have noted potential problems with the sky-subtraction 
in SDSS spectra at these wavelengths. Paris et al., for exam- 
ple, found a net positive residual below 4000 A in a sample 
of quasar spectra that should show no flux below 4280 A due 
to the presence of a DLA at z > 3.7. In principle, a positive 
residual would cause us to overestimate F(z) at z < 2.3, and 
would also impact our higher-redshift measurements, since 
the errors in our F(z) points will be correlated. We note 
that, for our measurements, the impact of sky-subtraction 
residuals in the far blue is likely to be small; we only use 
data at Aobs < 4000 A in the lowest-redshift composites 
(zcomp ^ 2.75), and these were constructed using relatively 
bright quasars. Nevertheless, we performed multiple tests to 
determine whether sky-subtraction errors could significantly 
affect our F(z) results. 

We first directly estimated the mean sky residual (and 
any contamination from intervening objects within the SDSS 
fibers) by measuring the flux over A b s = 3900 — 4000 A in 
the spectra of each of our quasars at z > 4.5. These objects 
should generally not have any flux below 4000 A due to 
the integrated Lyman-series opacity of the IGM and any 
optically thick absorbers. A few exceptionally clear lines 
of sight may show some transmission, in which case the 
flux will be an upper limit on the sky-subtraction resid- 
uals. We measured a mean flux of /4000 — (8 ± 1) x 
10~ 19 ergs -1 cm -2 A -1 . For our quasars at z < 2.8, this 
is, on average, roughly 0.5 per cent of the observed flux over 
the same wavelength range. This is comparable to (albeit 
slightly larger than) than the la errors on F(z)/F(z = 2.15) 
at z < 2.3, as described below, and substantially less than 
the errors in the final scaled F(z) values. This supports the 
conclusion that, since the quasars we use to measure the 
flux below 4000 A are bright, sky-subtraction errors at these 
wavelengths will have a relatively minor impact. 

Second, we tried adjusting the redshift interval over 
which we scaled our absolute F(z) values. If sky-subtraction 



errors are artificially increasing our F(z)/F(z = 2.15) val- 
ues at z < 2.3, the matching our results to literature values 
at these redshifts should yield a normalization factor that 
is too low. To test thi s , we t ried scaling our results to the 
IFaucher-Giguere et alj |2008h measurements over 2.6 ^ z ^ 
3.0, rather than over 2.2 ^ z ^ 2.5. The resulting scale fac- 
tor was higher by only 0.5 per cent, which is well within the 
la uncertainty of ~1 per cent (see measurements below). 

Finally, we tried fitting F(z) using only data at absorp- 
tion redshifts z > 2.5, or A D b s > 4250 A. This had a negli- 
gible impact on the shape of F(z) at z > 2.5. In summary, 
therefore, sky-subtraction residuals in the far blue do not 
appear to have a significant impact on our results. We em- 
phasize, however, that systematics such as these may need to 
be more carefully addressed in order to achieve substantially 
higher-precision mean flux measurements than the ones pre- 
sented here. 



3 RESULTS 

Our results for F( z)/F(z — 2.15), along with our F(z) 
values scaled to the IFaucher-Giguere et al.l (|2008h measure- 
ments at z ^ 2.5, are shown in Figure|4]and listed in Table(2] 
The errors are la and reflect only the diagonal elements of 
the covariance matrix. For F(z) the errors include the uncer- 
tainty in the overall absolute scaling, which has a fractional 
uncertainty of ~1 per cent (for the adopted set of corrections 
for continuum bias, metal-line contamination, and other sys- 
tematics made by Faucher-Giguere et al.) and dominates the 
error budget at z ^ 3.75. 

The mean flux exhibits a smooth decline with redshift 
over 2 < z < 5. Notably, we s ee no indication of th e fea- 
ture at z ~ 3.2 reported b y iBernardi et al.l (|2003T ) and 
IFaucher-Giguere et al.l (|2008l ). Comparisons with these and 
other literature values are discussed in Section [4] 

Our primary systematic uncertainty is whether the 
shape and amplitude of the mean unabsorbed continuum in 
the Lya forest is the same for all composites, regardless of 
redshift. To check this, we first reconstructed the continuum 
in the forest for each composite, scaled by F(z — 2.15), by 
dividing the observed fluxes by our raw (uncorrected for ab- 
sorption due to optically thick systems) relative transmitted 
flux values, 

/c 0nt (Ao, w P )^= 2 .i 5 ) = 4^;tTi5) ■ < 4 > 

Here, z is calculated from the observed wavelength of each 
pixel, and the transmitted flux is interpolated between the 
values plotted in Figure [4] Observed wavelengths for a com- 
posite at redshift z comp are related to the rest-frame wave- 
lengths plotted in Figure [2] simply by A b s = Ao(l + z C omp). 
We then computed the mean scaled continuum across all 
composites with z comp < 4.05, where the redshift limit was 
chosen to avoid noisier composites at higher redshifts. The 
result is shown in Figure [5] Finally, a model for the ob- 
served flux for each composite was created by multiplying 
the mean scaled continuum, / C ont(Ao)-F(.z = 2.15), by the 
relative mean transmitted flux, F(z)/F(z = 2.15), evalu- 
ated separately for each pixel. 

The fractional residuals from these models, (/ b s — 
/modei)//modoi, are shown in Figure [6] Error bars were com- 
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Figure 4. The mean transmitted flux in the Lyc? forest. Left panel: F(z) as a fraction of the mean transmitted flux at z = 2.15. These 
values are computed directly from our composite spectra, and include a correction for the abso r ption due to optically thick systems. 
Right panel: F(z) values after scaling our results to the measurements of iFaucher-Giguere et all l (200Sh over 2.2 ^ z ^ 2.5. Error bars 
are ltr, and include uncertainties in the scaling factor for F(z) points in the right-hand panel. Tabulated values are given in Table [2] 
The full covariance matrices can be found in the accompanying supplemental material. 
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Figure 5. The mean quasar continuum in the Lyc? forest. The 
continuum was calculated by averaging the composites with 
•Zcomp < 4.05 after correcting each one for the measured relative 
mean transmitted flux, F(z)/F(z = 2.15), calculated separately 
at each observed wavelength. The overall amplitude of the contin- 
uum has therefore implicitly been multiplied by F(z = 2.15). The 
solid line gives the results from the full sample. The shaded region 
shows the lc uncertainties calculated from bootstrap trials. 



puted by repeating the above procedure for each bootstrap 
trial. The absence of significant features in the residuals in- 
dicates that the unabsorbed continuum is highly similar be- 
tween composites. This includes the shape and amplitude of 
the emission lines as well as the overall slope of the contin- 
uum, which would introduce a tilt in the residuals if it varied 



significantly between composites. The self-similarity of the 
continuum supports our assumption that the ratios of the 
observed fluxes between composites can be modeled purely 
by the evolution in F(z) (plus a minor contribution from op- 
tically thick absorbers) , at least up to a single normalization 
factor. 

We emphasize that this test addresses the self-similarity 
of the continuum only over the Lycv forest. Variations with 
redshift in the observed mean quasar spectral energy distri- 
bution ( SEP) at other wavelengths ar e known to exist. For 
example. IWorseck fc Prochaskal |201ll ) found that SDSS, in 
order to avoid contamination from the stellar locus, prefer- 
entially selects quasars with red SEDs in the u and g bands 
over 2.7 < z < 3.5. They found that this primarily effects 
the observed mean flux for higher-order Lyman-series tran- 
sitions, and has relatively little impact on the Lya forest 
(see their Figure 17). Nevertheless, careful modeling of such 
selection effects may be required in order to achieve signifi- 
cantly more precise mean flux measurements than the ones 
presented here. 



3.1 Analytic function 

In addition to fitting a discreet F(z), we fit an analytic 
function to the effective optical depth as a function of 
redshift, Tes(z) = —lnF(z). It has been conventional to 
parametrize T e s as a simple power law of the form T e g(z) = 
to [(1 + z)/(l + zo)]' 8 . If one assumes this form for r e g(z), 
then both the normalization and the slope will be directly 
constrained from flux ratios, since neither to nor /3 cancels 
when computing F(zi)/ F(z2). With the precision and long 
redshift baseline provided by our dataset, however, we found 
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Figure 6. Fractional flux residuals in the Lyo forest for each of the twenty-six composites, plotted with respect to the corresponding 
Lya absorption redshifts. The residuals are calculated as (/ Q b s — /model)//modcl> where / mo dcl i s the product of the mean continuum 
shown in Figuro[5]and the relative mean transmitted flux, F(z)/F(z = 2.15), evaluated at each observed wavelength (i.e., redshift). ltr 
error bars were computed from bootstrap trials. A unity offset is added between composites for clarity. The lack of significant residuals 
demonstrates that the continuum over the Lya forest (i.e., the slope of the underlying power-law continuum and the relative shape and 
amplitude of the emission lines) changes little with quasar redshift, and that the ratios of observed fluxes between composites can be 
well modeled by the evolution in the mean Lya forest opacity. 
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Figure 7. Mean transmitted flux in the Lya forest expressed as an effective optical depth, r fr (z) = — ln_F(,z). The solid line is an analytic 
fit to the full data set (not the binned points displayed), using a modified power law of the form T e g- (z) = to [(1 + z)/(l + Zo)P +C. The 
shaded region gives the lc uncertainty in the model r c ff values, taking into account the covariance between the model parameters. See 
Section |3T] for details. The axes are scaled such that a pure power law with C = would appear as a straight line. 



that, although a pure power law could provide a reasonable 
fit to the observed flux ratios (i.e., produce an F(z) with the 
correct shape), the overall normalization of F(z) was sub- 
stantially too low when compared to literature values. We 
found, however, that a reasonable fit could be achieved by 
allowing for a constant offset to T c g(z), such that 

« W =^(i£L)\<7. (5) 

Following our approach for the discreet F(z), we fit to and /3 
directly using the flux ratios, then solved for C by scaling our 
results to those of IFaucher-Giguere et all (|2008l ) over 2.2 ^ 
z ^ 2.5. Bootstrapping was again used to compute error 
estimates, including the covariance between fit parameters. 

The results for T e g(z) are plotted in Figure]?] The best- 
fitting parameters are [r ,/3,C] = [0.751,2.90,-0.132], for 
zo = 3.5, and the full covariance matrix is given in Table [3] 
We emphasize that to and /3 were constrained directly from 
the composite flux ratios (corrected for absorption due to 
optically thick systems) rather than the binned F(z) values. 
Nevertheless, the analytic fit closely matches the discreet 
values at all redshifts once both are scaled to the literature 
values at 2 ~ 2. In F(z) space, the deviations of the discreet 
F(z) from the nominal analytic fit have a formal x 2 of 64, 
calculated using the full covariance matrix, for 28 redshift 
bins. This suggests that some of the residual structure in the 
discreet fit may be real, that a modified power law is not a 
sufficiently accurate model for T c g(z), or that we have some- 
what underestimated the errors in the discreet case. Even so, 
the departures of the discreet points from the analytic fit are 
small, with a mean fractional deviation in the flux of two per 
cent. A modified power law for T e a(z) thus appears to be a 



reasonable approximation to the true evolution of the Lya 
opacity over this redshift range, although we caution that 
the fit should not necessarily be extrapolated to higher or 
lower redshifts. 



4 COMPARISON TO PREVIOUS RESULTS 

We now provide a brief comparison to some previous re- 
sults from the literature. The most careful measurements 
of the mean Lya flux using high-resolution spectra has ar- 
guably been made by IFaucher-Giguere et ail (|2008l ). whose 
values at 2 2.5 we have used to normalize our own re- 
sults. Their values are plotted along with our results in Fig- 
ure [5] Although we normalized our values to theirs only 
over 2.2 ^ z ^ 2.5, the overall agreement out to z — 4.2 
is very good. We note that at z ^ 3.5 the Faucher-Giguere 
et al. values for T e g(z) are systematically somewhat higher 
(Ar c ff ~ 0.05) than ours. This may indicate that the contin- 
uum corrections made by those authors were somewhat too 
large at the high-redshift end. More significantly, the fea- 
ture in T e ff(z) at 2 ~ 3.2 tentatively identified by Faucher- 
Giguere et al. is absent from our results. In light of the 
substantially higher precision offered by the present analy- 
sis, the feature appears likely to have been spurious. It may 
either reflect random fluctuations in the Faucher-Giguere et 
al. data, or may be the result of a subtle systematic effect. 

A potential concern is whether an analysis using com- 
posite spectra, which average the flux from objects over a 
range of redshifts, may smooth out discreet features in the 
mean flux evolution. We tested this by generating artificial 
composites with the same redshift range and binning as used 
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Table 2. Tabulated values for the mean transmitted flux in the 
Lya forest. Columns give the redshift of each bin, the mean trans- 
mitted flux as a fraction of its value at z = 2.15 before and after 
correcting for absorption due to optically-thick absorbers, and the 
corrected abs olute mean transmitted flux af ter scaling to the mea- 
surements of iFau chcr-Gigu ere et all 1 120081 ) over 2.2 ^ z 2.5. 
The full covariance matrices can be found in the accompanying 
supplemental material. 
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Table 3. Covariance matrix for T e ff (z) parameters (equation [5} . 
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for the real data, i.e., bins of Az = 0.1 from z = 2.2 to 4.4, 
and Az = 0.4 from z — 4.4 to 5.2. Each composite was 
generated from 100 individual artificial spectra with red- 
shifts distributed uniformly across the bin. For simplicity, 
a uniform continuum was assumed for each spectrum. The 
input Lya forest abso rption for each individual sp ectrum 
was calculated from the lFaucher-Giguere et all l|200ct ) power 
law±Gaussian fit to r e a(z), which includes a Gaussian dip of 
width a z ~ 0.1 centere d at z ~ T his is comparable to 
the feature reported bv lBernardi et all (|2003T > . After adding 
noise to the artificial composites at a level comparable to 
that in the real data, we re-measured the mean flux evolu- 

See values using the ISchave et al.l ll2003h metal correction in 
Table 5 of lFaucher-Giguere et al.l d2008l~ 
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Figure 8. Our results as shown in F igure [7] plotted along with 
those of iFaucher-Giguere et a~ ] 11200811 , which arc the most care- 
ful measurements to date made with high-resolution spectra. Al- 
though our values were scaled to the Faucher-Giguere et al. mea- 
surements only at z ^ 2.5, the overall agreement out to z = 4.2 is 
very good. In light of the precision offered by the new measure- 
ments presented here, the feature at z ~ 3.2 tentatively identified 
by Faucher-Giguere et al. appears to be consistent with noise. We 
note that the Faucher-Giguere et al. points over 2.2 ^ z ^ 2.5, to 
which we normalized our results, appear noisy because they are 
presented on a log plot; their errors on the mean flux over this 
range are, in fact, relatively small (up < 0.02). 
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z 

Figure 9. A demonstration of the ability of our approach to re- 
co ver a feature in the T p ff (z) evo lutio n similar to the one rep orted 
bv lFaucher-Giguere etatl I I2008TI and lBernardi et al.l J2003 V The 
solid line represents the input r e fj (z), while the dotted line shows 
a pure power law for comparison. Points with error bars are the 
T e ff(z) we recover by applying out analysis technique to artifi- 
cial composite quasar spectra generated with the input T e ff(z) 
but with properties otherwise similar to those in our sample. See 
Section [4] for details. 

tion using the approach outlined in Section ^. 21 The results, 
plotted in Figure [9] show that the input flux evolution was 
recovered with high accuracy. Although this was a simpli- 
fied test, it demonstrates the key point that measurements 
of the mean flux evolution using composite spectra are en- 
tirely sensitive to discreet features similar to the ones that 
have been reported in the literature. 
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Figure 10. Our results as shown in Figure [7] plotted along with previous measurements made from SDSS spectra. Error bars have 
been omitted from the literature values for clarity and to highlight systematic trends. Deviations in the first three sets of measurements 
jBernardi et alj|2003l ; iMcDonald et al.ll2005l ; IDall'Aglio et alj|2009h from the results presented here arc likely to reflect systematic errors 
in those works. The overall agreement between the IPaTiset^i.1 ll201ll ) results and ours is generally good, although the Paris et al. 
measurements have not been corrected for metal-line contamination. See Section [4] for discussion. 



We next turn to previous measurements 
SDSS spect ra . Res ults from 
McDonald et all (120051 ) 



mg_ 



ma de us - 

Bernardi et alj (|2003h . 
IDall'Aglio et all < |2009h . and 



Paris et al.l ( 201 ll ) are plotted along with our measurements 
in Figure [lOl In the three earliest works there are systematic 
differences from our results. We propose that these differ- 
ences are most likely to be explained by systematic errors 
in the previous analyses, as discussed belo w. In contrast, 
the agreement between the lParis et al.l(|201ll ) measurements 
and our results is generally very good, which suggests that 
the measurements from SDSS spectra have now largely con- 
verged. 



iBernardi et al.l |200ot ) used composites generated from 
1061 individual quasar spectra to simultaneously constrain 
Tes(z) and the mean unabsorbed continuum over the Lya 
forest using two approaches. In the first approach, they 
conducted a joint fit to analytic functions for both T e g(z) 
and the mean continuum. Second, they used an iterative 
approach to fit a discreet T e g(z) along with an arbitrar- 
ily shaped continuum for each composite. Both approaches, 
however, face significant limitations. The composite spectra 
used by Bernardi et al., as well as the ones used here, fun- 
damentally do not contain enough information to constrain 



both the amplitude and shape of F(z) (or T e g (z)) unless ad- 
ditional assumptions are made. In their analytic approach, 
Bernardi et al. assumed a power la w + Gaussian form for 
T e g (z) similar to the one used by iFaucher-Giguere et al.l 
(120081). As discussed in Section 13.11 the ratios of fluxes at 



different redshifts directly provided by composite spectra 
will constrain both the normalization and slope of a pure 
power law. In this case, therefore, T e g(z) will be implicitly 
normalized purely by the choice of fitting function (an ad- 
ditional, narrow Gaussian component will not change this 
since the contribution to T e g(z) will go to zero far from the 
Gaussian centroid). Based on the results presented here, we 
would expect that the Bernardi et al. results would be off 
by at least a normalization factor, and indeed, an offset is 
seen in Figure \W\ that is consistent with the fact that we 
required a negative correction to T e fj (z) (C < 0) in order to 
correctly normalize our results. We note that there may be 
additional systematic errors related to the choice of analytic 
shape for the continuum, which Bernardi et al. parametrized 
as a power law plus three Gaussians to approximate the 
emission lines in the forest. 

For the secon d, iterative approach used by 
IBernardi et"afl (|2003T ), it is not clear that their results 
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were converged, as claimed. They used power laws nor- 
malized at 1450 A as initial estimates of the continua, 
and refined the continuum shapes based on the inferred 
F(z) calculated at each iteration. In tests with our data, 
we found that after 3-4 iterations, as used by Bernardi 
et al., this procedure roughly recovered the continuum 
shape plotted in Figure [5] However, the normalization of 
the continuum, and hence F(z), scaled directly with the 
normalization of the initial power-law continuum. Moreover, 
although the changes in the continuum were minor after 
a few iterations, the continua were not converged. After 
many iterations, the continua evolved to incorporate all of 
the relative differences between composites, producing a 
constant F(z). This appears to be a generic result when 
the continua are allowed to converge separately for each 
composite. If the continua were held to be the same for 
each composite, then the shape of the continuum did 
converge, producing a relative F(z) very similar to the 
one presented here. The overall normalization, however, 
did not converge. This simply demonstrates the fact that, 
without a completely unabsorbed composite to serve as 
a reference, the composites themselves do not contain 
enough information to fully constrain both F(z) and the 
unabsorbed continuum without making other assumptions. 

iDaU'Aglio et all (|2009l ) fitted continua to 1733 individ- 
ual quasar spectra by combining localized spline fits with 
a redshift-dependent global correction estimated from mock 
sp ectra. This approach is i n prin ciple similar to the one used 
bv lFaucher-Giguere et all (|2008T ); however, there are at least 
two significant differences. First, the use of low-resolution 
data means that the corrections must be substantially larger, 
and will be significant even down to z ~ 2. Additionally, the 
mock spectra used by DalPAglio et al. were created by ran- 
domly generating absorption lines from fixed number den- 
sity, column density, and Doppler parameter distributions, 
rather than produced from cosmological simulations, as used 
by Faucher-Giguere et al.. The continuum corrections esti- 
mated from these mock spectra will therefore be subject to 
inaccuracies in the line parameter distributions, which were 
themselves estimated from high-resolution spectra that may 
have suffered from continuum fitting errors. This systematic 
effect offers a reasonable explanation for the relative tilt in 
the evolution of r e s(z) measured by DalPAglio et al. with 
our respect to our values. We note, however, that DalPAglio 
et al. find a smooth evolution with redshift, similar to our 
results, w ithout any evid e nce fo r the feature at z ~ 3.2 re- 
por ted byfeernardi et all (|2003h . 

iMcDonald et all (120051 ) also found a relatively smooth 
evolution of T e ff (z) with redshift in a sample of 3035 quasars. 
We include this work for completeness, although the authors 
stress that their method for estimating individual quasar 
continua in the Lya forest using principle component analy- 
sis was not necessarily converged. They also required a nor- 
malization factor taken from other measurements, though 
this would not affect the overall shape of r c f((z). In gen- 
eral, their results reinforce the general conclusion that the 
mean transmitted flux measured from SDSS spectra evolves 
smoothly with redshift, consistent with the results presented 
here. 

Finally, iParis et all (120111 ) recently applied PCA to 
measure the mean transmitted Lya flux in the SDSS spectra 
of 2576 quasars. Their approach used principle components 



derived from z ~ 3 quasar spectra to predict the unab- 
sorbed continuum in the Lya forest of individual objects 
based on the flux redward of Lya. Overall, their values are 
generally consistent with our results, although very mild sys- 
tematic differences can be seen in Figure 1101 We caution 
that, despite the apparent agreement, a direct comparison 
with their results may not be entirely appropriate. Although 
Paris et al. addressed the issue of optically thick systems by 
avoiding lines of sight containing DLAs, they did not cor- 
rect for metal line contamination. Since they were measuring 
the mean flux directly from continuum-normalized spectra, 
such a correction would have been necessary to calculate 
the mean transmitted flux due to Lya alone. We note that 
applying an adjustment of Ar c ft ~ —0.03 for metals (see 
Section [2. 2. 1|) would bring their t c h(z) values systematically 
below our one-sigma limits everywhere except in their three 
highest redshift bins. 



5 SUMMARY 

We have presented new measurements of the mean trans- 
mitted flux in the Lya forest over 2.15 ^ z ^ 4.85 made 
using composite quasar spectra. Our sample includes 6065 
individual spectra from the SDSS DR7 quasar catalogue, 
which were combined into 26 composites with mean quasar 
redshifts z comp — 2.25 to 5.08. 

The fact that the quasar composites all have highly sim- 
ilar fluxes redward of Lya suggests that the differences in 
the observed flux in the forest are driven by the evolution in 
the mean intergalactic Lya opacity. The ratio of F(z) at two 
different redshifts can therefore be measured by comparing 
the observed mean flux in two composites a the same rest 
wavelength. We used these ratios to directly constrain F(z) 
as a fraction of F(z = 2.15). The overall normalization for 
F(z) was then determined using literature values made from 
high-resolution data at z ^ 2.5, where continuum uncertain- 
ties, although accounted for, are relatively minor. The sta- 
tistical accuracy offered by the quasar spectra means that 
this normalization factor now dominates the uncertainty in 
F(z) at z < 3.75. 

Our main results give the mean transmitted flux associ- 
ated with systems of column density JVhi < 10 19 cm -2 , and 
do not include the contribution from metal lines. We fit the 
mean flux using two functions. The first is a discreet F(z) 
in bins of Az = 0.1, while the second is a modified power 
law of the form T e s(z) = to [(1 + z)/(l + zo)] 13 + C. where 
r e ff(z) = — \nF(z). Although previous works had generally 
parametrized r e ff(z) as a pure power law (C = 0), we found 
that this was insufficient to describe the observed flux ratios 
at the required level of precision while simultaneously pro- 
ducing the correct normalization. We emphasize that the 
true functional form of t c s{z) is unknown, and that our 
parametrization is appropriate only in that it provides a rea- 
sonable match to the discreet F(z) values. Caution should 
be used when extrapolating this (or any) fit to higher and 
lower redshifts. We also stress that, due to the fact that 
our results are generated from the ratios of fluxes at differ- 
ent redshifts, the parameters for both the discreet F(z) and 
modified power law T e g{z) are correlated, and that the full 
correlation matrices should be considered when applying our 
results. 
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The mean transmitted flux is found to evolve smoothly 
with redshift, which is consistent with a gradual evolu- 
tion of the H I ionization rate and gas temperature over 
2 < z < 5. We see no evidence of previously reported 
features at z ~ 3.2, which appear to have been spurious. 
Our results are otherwise generally consistent with the most 
caref ul measurements made to d ate using high-resolution 
data |Faucher-Giguere et al.ll2008j ) , but offer a substantially 
greater level of precision and extend out to higher redshifts. 
Our measurements also offer an improvement over previous 
efforts using SDSS spectra, which appear gener ally to have 
been affected by systematic errors (although see lParis et al.l 

Soil). 

This work demonstrates the significant advantages of 
using large data sets to measure the mean transmitted flux. 
In addition to the high degree of statistical accuracy, the use 
of composites makes it possible to avoid fitting continua at 
high redshifts, which is the primary source of systematic un- 
certainty in other approaches. These measurements should 
therefore serve as a reference for a variety of studies that 
make use of the Lya forest. 
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